R version 4.0.3 (2020-10-10) – “Bunny-Wunnies Freak Out”
Packages used for NMDS: vegan (version 2.5-7)
The document shows a series of NMDS ordinations for reference benthic communities in Virginia with environmental characteristics overlaid to evaluate natural differences in community compositions across Virginia. These NMDS will support the Genus level IBI development process. This analysis is the first run of all of reference sites and sites that were deemed reference piedmont sites in Virginia. No West Virginia DEP data is used in this analysis.Reference sites were evaluated by regional biologists.
The dataset used includes all reference stations collected in Virginia that were deemed reference through a series or water quality parameter filters and regional biologist review. If stations appeared in the dataset more than 4 times, then the most recent 4 samples were used and the rest removed. Samples that had a total number of taxa below 100 collected at the time of sampling were also removed. Taxa that occurred in the dataset <= 5% of the time were removed. The data was log10 +1 transformed. Environmental factors were compiled for each station and used to plot over the NMDS to show environmental variation associated with the community matrix. The envfit function in Vegan was used to plot the continuous environmental variables. Some environmental variables like precipitation, slope, and elevation have not been calculated for all watersheds yet and will be added at a later date.
The first step was to read in the reference site bug taxa list and environmental factors dataset for each station. Join the environmental dataset with the bug dataset to account for multiple observations of each station and collection date and time.
Check to make sure the bug and environmental join was successful:
Number of rows in Community Matrix: 783
Number or rows in Environmental Matrix: 785
The data was log10+1 transformed. Rare taxa (<=5%) were removed.
## Run 0 stress 0.1787787
## Run 1 stress 0.1790783
## ... Procrustes: rmse 0.002504029 max resid 0.05042826
## Run 2 stress 0.1800402
## Run 3 stress 0.1791169
## ... Procrustes: rmse 0.002979715 max resid 0.0506303
## Run 4 stress 0.1790946
## ... Procrustes: rmse 0.003184204 max resid 0.04584174
## Run 5 stress 0.1796415
## Run 6 stress 0.1794197
## Run 7 stress 0.1790037
## ... Procrustes: rmse 0.002487529 max resid 0.04534014
## Run 8 stress 0.179314
## Run 9 stress 0.178976
## ... Procrustes: rmse 0.002563227 max resid 0.04113389
## Run 10 stress 0.1789059
## ... Procrustes: rmse 0.002089364 max resid 0.04161214
## Run 11 stress 0.1870266
## Run 12 stress 0.1792536
## ... Procrustes: rmse 0.003577928 max resid 0.05002788
## Run 13 stress 0.1788356
## ... Procrustes: rmse 0.001888102 max resid 0.04268624
## Run 14 stress 0.1846527
## Run 15 stress 0.180944
## Run 16 stress 0.1789466
## ... Procrustes: rmse 0.001634449 max resid 0.04524993
## Run 17 stress 0.1809055
## Run 18 stress 0.1790442
## ... Procrustes: rmse 0.00266505 max resid 0.04579269
## Run 19 stress 0.1793592
## Run 20 stress 0.1790002
## ... Procrustes: rmse 0.002942282 max resid 0.05073364
## Run 21 stress 0.1866547
## Run 22 stress 0.1790225
## ... Procrustes: rmse 0.00298488 max resid 0.04545408
## Run 23 stress 0.1791671
## ... Procrustes: rmse 0.003386359 max resid 0.04587383
## Run 24 stress 0.1788438
## ... Procrustes: rmse 0.001474478 max resid 0.03801694
## Run 25 stress 0.1790461
## ... Procrustes: rmse 0.002782981 max resid 0.05032982
## Run 26 stress 0.1790258
## ... Procrustes: rmse 0.003040645 max resid 0.04547114
## Run 27 stress 0.1790523
## ... Procrustes: rmse 0.002707138 max resid 0.04568219
## Run 28 stress 0.1790441
## ... Procrustes: rmse 0.002661951 max resid 0.04579274
## Run 29 stress 0.1793933
## Run 30 stress 0.1790974
## ... Procrustes: rmse 0.003154143 max resid 0.04549251
## Run 31 stress 0.1810013
## Run 32 stress 0.1790062
## ... Procrustes: rmse 0.002335555 max resid 0.04532854
## Run 33 stress 0.1789101
## ... Procrustes: rmse 0.001911388 max resid 0.05045323
## Run 34 stress 0.1788858
## ... Procrustes: rmse 0.001930166 max resid 0.03760034
## Run 35 stress 0.1789635
## ... Procrustes: rmse 0.002407358 max resid 0.05052184
## Run 36 stress 0.1789768
## ... Procrustes: rmse 0.002607093 max resid 0.04162671
## Run 37 stress 0.1790041
## ... Procrustes: rmse 0.003014302 max resid 0.05080422
## Run 38 stress 0.1791678
## ... Procrustes: rmse 0.002934758 max resid 0.0500105
## Run 39 stress 0.1789644
## ... Procrustes: rmse 0.002435676 max resid 0.05055896
## Run 40 stress 0.1792212
## ... Procrustes: rmse 0.003240301 max resid 0.05010465
## Run 41 stress 0.1791985
## ... Procrustes: rmse 0.003305406 max resid 0.0498025
## Run 42 stress 0.1788399
## ... Procrustes: rmse 0.001444341 max resid 0.03847479
## Run 43 stress 0.1797901
## Run 44 stress 0.1810849
## Run 45 stress 0.1788975
## ... Procrustes: rmse 0.002097443 max resid 0.04075255
## Run 46 stress 0.1789609
## ... Procrustes: rmse 0.002398921 max resid 0.05052494
## Run 47 stress 0.1788167
## ... Procrustes: rmse 0.00153419 max resid 0.0345005
## Run 48 stress 0.179167
## ... Procrustes: rmse 0.002921762 max resid 0.05000944
## Run 49 stress 0.1797454
## Run 50 stress 0.1788138
## ... Procrustes: rmse 0.001338495 max resid 0.03340206
## Run 51 stress 0.1790042
## ... Procrustes: rmse 0.002157104 max resid 0.0452859
## Run 52 stress 0.1789642
## ... Procrustes: rmse 0.002428818 max resid 0.05054606
## Run 53 stress 0.1790493
## ... Procrustes: rmse 0.00256716 max resid 0.04542099
## Run 54 stress 0.180569
## Run 55 stress 0.1790144
## ... Procrustes: rmse 0.002950913 max resid 0.04240262
## Run 56 stress 0.179546
## Run 57 stress 0.1790032
## ... Procrustes: rmse 0.002924426 max resid 0.05071931
## Run 58 stress 0.1799167
## Run 59 stress 0.1788972
## ... Procrustes: rmse 0.002078527 max resid 0.04070022
## Run 60 stress 0.1793426
## Run 61 stress 0.1794974
## Run 62 stress 0.1868693
## Run 63 stress 0.1793748
## Run 64 stress 0.1793932
## Run 65 stress 0.178984
## ... Procrustes: rmse 0.00232877 max resid 0.0453998
## Run 66 stress 0.1806733
## Run 67 stress 0.1799345
## Run 68 stress 0.1788164
## ... Procrustes: rmse 0.001609942 max resid 0.03542758
## Run 69 stress 0.1794745
## Run 70 stress 0.1789966
## ... Procrustes: rmse 0.002265871 max resid 0.04567002
## Run 71 stress 0.1802257
## Run 72 stress 0.1790471
## ... Procrustes: rmse 0.002737429 max resid 0.05029246
## Run 73 stress 0.1802009
## Run 74 stress 0.1789141
## ... Procrustes: rmse 0.001894929 max resid 0.05046208
## Run 75 stress 0.1789612
## ... Procrustes: rmse 0.002405158 max resid 0.05053038
## Run 76 stress 0.1790018
## ... Procrustes: rmse 0.002169604 max resid 0.04532874
## Run 77 stress 0.1791345
## ... Procrustes: rmse 0.003035987 max resid 0.0506646
## Run 78 stress 0.1789539
## ... Procrustes: rmse 0.002515393 max resid 0.05065274
## Run 79 stress 0.1833642
## Run 80 stress 0.1790959
## ... Procrustes: rmse 0.003114367 max resid 0.05006005
## Run 81 stress 0.1799673
## Run 82 stress 0.1788391
## ... Procrustes: rmse 0.00157813 max resid 0.04073087
## Run 83 stress 0.1790027
## ... Procrustes: rmse 0.002924881 max resid 0.05072921
## Run 84 stress 0.1800514
## Run 85 stress 0.1789921
## ... Procrustes: rmse 0.002353372 max resid 0.05019188
## Run 86 stress 0.1790622
## ... Procrustes: rmse 0.003385237 max resid 0.05096852
## Run 87 stress 0.1794752
## Run 88 stress 0.1792028
## ... Procrustes: rmse 0.00332282 max resid 0.04974752
## Run 89 stress 0.178955
## ... Procrustes: rmse 0.00252075 max resid 0.05067149
## Run 90 stress 0.1789136
## ... Procrustes: rmse 0.001882456 max resid 0.05044315
## Run 91 stress 0.1800532
## Run 92 stress 0.1799545
## Run 93 stress 0.1790093
## ... Procrustes: rmse 0.002213583 max resid 0.04533748
## Run 94 stress 0.1797498
## Run 95 stress 0.1790201
## ... Procrustes: rmse 0.002997493 max resid 0.04546739
## Run 96 stress 0.1791602
## ... Procrustes: rmse 0.003373377 max resid 0.05073323
## Run 97 stress 0.1797267
## Run 98 stress 0.1790374
## ... Procrustes: rmse 0.003197673 max resid 0.0509233
## Run 99 stress 0.1800704
## Run 100 stress 0.1790026
## ... Procrustes: rmse 0.002303792 max resid 0.04573901
## Run 101 stress 0.1789469
## ... Procrustes: rmse 0.001639198 max resid 0.04524654
## Run 102 stress 0.1793303
## Run 103 stress 0.1791304
## ... Procrustes: rmse 0.002928179 max resid 0.05054952
## Run 104 stress 0.1791309
## ... Procrustes: rmse 0.003002259 max resid 0.05063083
## Run 105 stress 0.179005
## ... Procrustes: rmse 0.00218341 max resid 0.04529602
## Run 106 stress 0.1791979
## ... Procrustes: rmse 0.003283326 max resid 0.04983816
## Run 107 stress 0.1806952
## Run 108 stress 0.1788983
## ... Procrustes: rmse 0.002111765 max resid 0.04074243
## Run 109 stress 0.179694
## Run 110 stress 0.1810268
## Run 111 stress 0.1791611
## ... Procrustes: rmse 0.00345752 max resid 0.05080635
## Run 112 stress 0.1834986
## Run 113 stress 0.1792974
## Run 114 stress 0.1791131
## ... Procrustes: rmse 0.003026405 max resid 0.05066145
## Run 115 stress 0.1789099
## ... Procrustes: rmse 0.001898634 max resid 0.05043461
## Run 116 stress 0.1789104
## ... Procrustes: rmse 0.001921786 max resid 0.05046668
## Run 117 stress 0.1795866
## Run 118 stress 0.1820863
## Run 119 stress 0.1799267
## Run 120 stress 0.1794766
## Run 121 stress 0.1788353
## ... Procrustes: rmse 0.001882728 max resid 0.04270471
## Run 122 stress 0.1793657
## Run 123 stress 0.1790809
## ... Procrustes: rmse 0.002485292 max resid 0.05044364
## Run 124 stress 0.1793388
## Run 125 stress 0.1789624
## ... Procrustes: rmse 0.002379664 max resid 0.05051404
## Run 126 stress 0.1818149
## Run 127 stress 0.179439
## Run 128 stress 0.178884
## ... Procrustes: rmse 0.001988346 max resid 0.0404745
## Run 129 stress 0.1815045
## Run 130 stress 0.1796759
## Run 131 stress 0.1797622
## Run 132 stress 0.1788563
## ... Procrustes: rmse 0.002466209 max resid 0.04311143
## Run 133 stress 0.1788858
## ... Procrustes: rmse 0.001962609 max resid 0.03725129
## Run 134 stress 0.1818615
## Run 135 stress 0.17944
## Run 136 stress 0.181253
## Run 137 stress 0.1789997
## ... Procrustes: rmse 0.00292363 max resid 0.05071451
## Run 138 stress 0.1795222
## Run 139 stress 0.1800116
## Run 140 stress 0.1792009
## ... Procrustes: rmse 0.003273409 max resid 0.04981916
## Run 141 stress 0.1871008
## Run 142 stress 0.179262
## ... Procrustes: rmse 0.003228282 max resid 0.05091956
## Run 143 stress 0.1805659
## Run 144 stress 0.1791994
## ... Procrustes: rmse 0.003329148 max resid 0.04976324
## Run 145 stress 0.1795558
## Run 146 stress 0.186486
## Run 147 stress 0.1797707
## Run 148 stress 0.1855824
## Run 149 stress 0.178953
## ... Procrustes: rmse 0.002463701 max resid 0.05059836
## Run 150 stress 0.1850121
## Run 151 stress 0.1789631
## ... Procrustes: rmse 0.002379757 max resid 0.05049656
## Run 152 stress 0.179668
## Run 153 stress 0.1788292
## ... Procrustes: rmse 0.001874587 max resid 0.03740593
## Run 154 stress 0.1791609
## ... Procrustes: rmse 0.003380909 max resid 0.05073613
## Run 155 stress 0.1861242
## Run 156 stress 0.1794818
## Run 157 stress 0.1797745
## Run 158 stress 0.1791128
## ... Procrustes: rmse 0.003014708 max resid 0.05064558
## Run 159 stress 0.1792149
## ... Procrustes: rmse 0.003888058 max resid 0.05109601
## Run 160 stress 0.1790411
## ... Procrustes: rmse 0.002625705 max resid 0.04578362
## Run 161 stress 0.179157
## ... Procrustes: rmse 0.003361845 max resid 0.05071439
## Run 162 stress 0.1787763
## ... New best solution
## ... Procrustes: rmse 0.0003581806 max resid 0.0085787
## ... Similar to previous best
## *** Solution reached
##
## Call:
## metaMDS(comm = NoncoastalFive[, 6:108], k = 3, trymax = 1000)
##
## global Multidimensional Scaling using monoMDS
##
## Data: NoncoastalFive[, 6:108]
## Distance: bray
##
## Dimensions: 3
## Stress: 0.1787763
## Stress type 1, weak ties
## Two convergent solutions found after 162 tries
## Scaling: centring, PC rotation, halfchange scaling
## Species: expanded scores based on 'NoncoastalFive[, 6:108]'
## NMDS1 NMDS2 r2 Pr(>r)
## Year 0.99272 0.12042 0.0160 0.21
## JulianDate -0.40251 0.91542 0.7186 0.01 **
## Latitude -0.54298 -0.83975 0.0221 0.10 .
## Longitude -0.63370 -0.77358 0.0267 0.06 .
## totalArea_sqMile -0.77223 -0.63534 0.3317 0.01 **
## ELEVMEAN 0.92445 0.38130 0.2053 0.01 **
## SLPMEAN 0.89201 0.45201 0.1469 0.01 **
## wshdRain_mmyr -0.76226 -0.64727 0.3131 0.01 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Permutation: free
## Number of permutations: 99
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Season, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Fall Spring
## delta 0.6362 0.6316
## n 396 387
##
## Chance corrected within-group agreement A: 0.046
## Based on observed delta 0.6339 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$X, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Ref Ref-Pied
## delta 0.6547 0.6492
## n 684 99
##
## Chance corrected within-group agreement A: 0.01578
## Based on observed delta 0.654 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$US_L3NAME, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Blue Ridge Central Appalachians Northern Piedmont Piedmont
## delta 0.5946 0.6055 0.6492 0.6329
## n 154 40 124 213
## Ridge and Valley
## delta 0.6573
## n 252
##
## Chance corrected within-group agreement A: 0.04529
## Based on observed delta 0.6344 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Basin_Code, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Appomattox Chowan-Dismal James-Middle James-Upper New Potomac-Lower
## delta 0.5805 0.584 0.6296 0.6532 0.6506 0.6493
## n 18 24 81 114 78 40
## Potomac-Shenandoah Rappahannock Roanoke Tennessee-Big Sandy
## delta 0.6553 0.6397 0.62 0.6036
## n 41 128 101 16
## Tennessee-Clinch Tennessee-Holston Yadkin York
## delta 0.6268 0.6257 0.5709 0.6224
## n 50 53 4 35
##
## Chance corrected within-group agreement A: 0.04623
## Based on observed delta 0.6338 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$ASSESS_REG, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## BRRO NRO PRO SWRO VRO
## delta 0.646 0.6644 0.619 0.6334 0.6481
## n 258 197 55 156 117
##
## Chance corrected within-group agreement A: 0.02704
## Based on observed delta 0.6465 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##Bioregion: Non-Coastal
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Bioregion, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Mountain Piedmont
## delta 0.6473 0.6493
## n 446 337
##
## Chance corrected within-group agreement A: 0.02454
## Based on observed delta 0.6482 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Gradient, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## MACS Riffle
## delta 0.6094 0.6587
## n 39 744
##
## Chance corrected within-group agreement A: 0.01239
## Based on observed delta 0.6563 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Order, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## 1 2 3 4 5
## delta 0.6367 0.6589 0.6305 0.6119 0.5699
## n 266 216 165 103 33
##
## Chance corrected within-group agreement A: 0.04378
## Based on observed delta 0.6354 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$StreamCate, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Large Medium Small
## delta 0.6043 0.6305 0.6536
## n 136 165 482
##
## Chance corrected within-group agreement A: 0.03665
## Based on observed delta 0.6401 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$WQS_CLASS, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## III IV V VI
## delta 0.6513 0.6448 0.6313 0.6012
## n 309 184 87 203
##
## Chance corrected within-group agreement A: 0.04507
## Based on observed delta 0.6346 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$WQS_TROUT, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## Yes
## delta 0.6598 0.6221
## n 493 290
##
## Chance corrected within-group agreement A: 0.02808
## Based on observed delta 0.6458 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
#Ecoregion separated by Season
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$BioregionSeason, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## MountainFall MountainSpring PiedmontFall PiedmontSpring
## delta 0.6234 0.6052 0.6091 0.6241
## n 225 221 171 166
##
## Chance corrected within-group agreement A: 0.07403
## Based on observed delta 0.6153 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
##
## Call:
## mrpp(dat = bugsnms_noncoast[, 6:108], grouping = samplescoresenv_noncoast$Bioregionsize, distance = "bray")
##
## Dissimilarity index: bray
## Weights for groups: n
##
## Class means and counts:
##
## LargeMountainFall LargeMountainSpring LargePiedmontFall
## delta 0.4825 0.5487 0.573
## n 29 27 42
## LargePiedmontSpring MediumMountainFall MediumMountainSpring
## delta 0.588 0.5731 0.5653
## n 38 51 50
## MediumPiedmontFall MediumPiedmontSpring SmallMountainFall
## delta 0.5818 0.5952 0.6105
## n 32 32 145
## SmallMountainSpring SmallPiedmontFall SmallPiedmontSpring
## delta 0.5768 0.5983 0.6072
## n 144 97 96
##
## Chance corrected within-group agreement A: 0.1192
## Based on observed delta 0.5853 and expected delta 0.6645
##
## Significance of delta: 0.001
## Permutation: free
## Number of permutations: 999
#Bioregion separated by Season